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I. 



INTRODUCTION 



The Commander in Chief Pacific has the basic responsi- 
bility of exercising control of the military forces of the 
United States from the west coast of North America to the 
east coasts of Asia and Africa. Since most of this vast 
area is either ocean or accessible by sea lines of communi- 
cation, the U.S. Navy has a critical role in the maintenance 
of a U.S. military presence in the Pacific Command area. 

The problem addressed in this thesis was posed by Vice 
Admiral E.P. Aurand, USM to Karl Eulenstein^ Senior Analyst in 
Section J 77 at CINCPAC. As COMTKIRDFLT, Admiral Aurand 
wanted to know how many aircraft carriers would be required 
to patrol the western littoral of the Pacific Ocean Basen 
as well as the littorals of the Indian Ocean Basen in order 
that at least one patrolling carrier would be within a 
command specified distance of a randomly occurring trouble 
spot with at least a command specified probability. More 
formally, what is the probability distribution of the minimum 
of the distances between n carriers and a single random 
trouble spot somewhere along the patrolled coast 

This thesis addresses the analytical expression of, 

P^[D = min(D^,D2, . . . ,D^) £d] (1) 
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where : 



n is the number of carriers on station 

D, ,D« ... D are the distances from each carrier to 
12 n 

a random trouble spot 

d is a command specified distance measured from, 
the trouble spot. 

The model allows the determination of the minimum n required 
in expression (1) such that the probability is at least as 
large as p. Other related inferences and applications are 
discussed. 

In the following a one-dimensional model and two approxi- 
mations that give useful bounds to it are developed. Then 
the analogous two-dimensional bounding models are considered. 
Finally the concept of a "sphere of influence" is incorporated 
into the one-dimensional model. 
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II. INITIAL ANALYTIC SOLUTION 



A. ASSUMPTIONS 

To increase the tractability of the analysis several 
simplifying assumptions were made: 

(1) The complex curvilinear coastlines v/hich mark the 
western boundry of CINCPAC'S area of responsibility are 
approximated by a straightline of length L. 

(2) Each of the n carriers were assigned individual 
patrol line segments of length ^ miles. 

(3) A carrier's position is assumed to be uniformly 
distributed on its assigned patrol line segment. 

(4) A trouble spot will appear at random along the 
coast line. 

(5) An individual carrier is considered as a point. 

(The effects the carrier airwing's radius of influence are 
analysed in Chapter V.) 

(6) Each carrier's location within its patrol segment 
is independent of the locations of the other n-1 carriers . 
The location of a trouble spot T is independent of the loca- 
tions of the carriers. See Figure 1. 

B . DEVELOPMENT 

The Mean Value Theorem for Integrals is utilized in the 
derivation of the analytical results detailed below. This 
theorem shows one may approximate a probability density 
function at a point h as 



12 



Schematic and Symbology for the Analytic Model 
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L 



X / i ~ 1 / 2 / 
D^, i = 1,2, 



L 

1 

n 



n the location of carrier i within 
line segment i 

n the distance to carrier i measured 
from a randomly occurring trouble 
spot t 

length of line to be patrolled 

the number of carriers assigned to 
patrol a line of length L 



t 



randomly occurring trouble spot 



^Notationally this paper will follow the custom of using 
uppercase letters for random variables and lowercase letters 
for parameters and outcomes with the exception of the 
parameter L which will denote the length of the line to 
be patrolled. 



Figure 1 



f (h) Ah 



= lim P[Xe (h,h+Ah) ] . 

Ah->0 + 

Examination of the geometry of the problem suggests 
that the development of the c\amulative distribution function, 
Fj^(d)^ = P^[D = min , . . . , D^) £d] be segmented as follows 

(1) Left end segment such that T e ^ ^ 

(2) Middle segments such that T e ( for some 

i = 2 , 3 , . - . , n-1 ^ ^ 

(3) Right end segment such that T e ^^,L) 
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Figure 2. Line segmentation 

The general analytic result was developed by conditioning 
on the segment containing the target. As a result of the 
assumption that the trouble spot will appear at random along 
the line (0,L), the target is located in segment i with 



The cumulative distribution function which will be 
developed is actually a function of the minimum of the 
distances, d; the length of the patrolled line, L; and the 
number of patrolling units, n; Fj^ (d; L,n) . For simplicity 
this terminology will be abbreviated as Fj^(d). 
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probability 1/n; that is 




for i = 1,2 



• • • / n • 



( 2 ) 



Given that T e [ 



(i-l)n in 
L ' L 



— ] , we further condition on T being 

Ij 



located in the first or second half of this line segment. 



D = [min(D^_^,D^ °i+l^ ^ ^ random variable 



whose cdf is being developed. The conditional pdf of D 
given T is unconditional over T (being located in a specific 
half of a line segment) and integrated over the range (0,d) 
to determine the conditioned cumulative distribution function 
(cdf) of D = [min ,D^_j_^) ] given that T is an element 

of a specific line segment i. 

1 . Left-hand Segment 

The assumptions made in developing this portion of 
the model are: 

(1) T '"Uniform (0,L), so Ts (0,— ) with probability 



1/n 

(2) ~ Uniform (0/^) 

(3) X„ Uniform ( , ill ) 

^ n n 

(4) X 2 / T are independent random variables, 
a. Conditional Probability Density Function 

Let S(i) be the event that the random variable 




n 



th 

T occurs in the i segment ( 



(i-1) L iL 
n ' n 



— ) of the line; then 



n 



P[S(i) ] = 1/n. 
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Consider the case t e 

D = = 1T-X^|, 

Zj 

and the maximum value of D is ^ ~ 0 £ d £ t, 

^D|T,S(i) = P [D £ (d,d+Ad) ] 

= P[X^ £ ( t-d, t-d-Ad) or X^ £ (t+d, t+d+Ad) ] 

= 2Adn/L . 

For t _< d £ ^ - t, 

P [D £ (d,d+Ad) ] = P [X^ £ (t+d, t+d+Ad) ] = 

Consider next the case t £ (^f~) , so D = min ( 0 ^, 02 ) 
D £ t . For 0 £ d £ ^ - t , 

P [D £ (d,d+Ad)] = P [X^ £ (t-d-Ad, t-d) 

or Xj^ £ (t+d, t+d+Ad) ] = 2Adn/L . 



Adn/L . 



and 
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For 



H 1 <1 i 



t t 



P[D e (d,d+Ad) ] = P[{X^ e (t-d-Ad,t-d) and > t+d+Ad} 

or {X^ < t-d-Ad and X^ £ (t+d,t+d+Ad) } ] 



= (Adn/L) • - t - d) . J + ( (t-d) -g) • (^) 

= Ad ^ t - d - Ad + (t-d-Ad) ] 

It n 

= 2Adn/L • [1 - . 

Thus it is seen that; 



folT^lt) 



2n/L ; 0 i <i < t 

n/L ; t £ d £ (^) - t 

f,- 0<d£(t)-t 

f(l-f): (|)-t<d<t 



for t £ (0,^) 

for t£(^,^) 
(3) 
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b. Unconditioning Over te 
t e (^,^) Given S(l) 

The conditional probability density function is 
graphed below, where values of fj^j^(dlt) are shown for various 
regions of the d-t plane - 




Figure 3. Unconditional probability density function 
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We now find the conditional distribution of D, 
given S(.l), by taking the expectation of f j-, | ,p S(l) 
with respect to the conditional distribution of T given 
S(l). Thus unconditioning over T (given S(l)) has the 
mathematical form: 

L/n 

^d|s(.1) ^ ^D|T,S(1) ^T|S(1)^^^ • 

0 



For 0 < d < L/2n, 



^D| S(l) 



d 

/ r 






gdt. 



2n 

L 






L 

n 






H-'l 






+ 2(1)^(£- 2d) + 2(f)^(l -^) (i 

Li ij n Li Li n 



|td) 



.n,, 2 ,2L , _ ,2 , , 

= (t-) [— d - 2nd /L] 

■Ij n 



^ < <3 1 J - 






L 

n 



^DlS(l) 



n n , r " 2n,, nd. n ,, 



= (^)2 (t-d)+2(g)^(l-^^) (L-d) 



= (r) ^ (— - 5d + 2d^n/L) . 

Li n 
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Therefore in the lefthand segment, 



^ - d - 2nd^/L] ; for 0 < d < ^ 

ii n — — . 

(?) ^ - 5d + 2d^n/L] ; for :^ < d < 

ij n zn — — 

(4) 

c. The Conditional Cumulative Distribution Function 
The cumulative distribution function (cdf) of 
D conditioned by S(l) in accordance with expression (2) is 
calculated by integrating the pdf over the range (0,d) of 
D: 



■dIs (1) 



(d) = 



d 

0 



For 0 < d < T?- 
— — 2n 



/ 



d 

^d|s( 1)^^^ ~ ^ [— -h-2nh /L] dh 

0 



,n^2f2Ld _ d^ _ 2nd^i 
^ n 2 3L ^ 



and 



F 



D I T^2n^ 



ii 

24 
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For ^ < d < ^ . 
Zn — — n 



Fr^lm(d t) 

D T ' 



L/2n 



/ (J)'l^-5h.2h2gl ah . 1 | 



,n^2 r3Ld 5d^ 2nd^ , ,5^,L^2 

^ n " 2 3L " ^ n ^ ^2^ 



- (2n\ /^\ 3 ■ 2 ^ 

^3L^ ^2n^ 24 



,n,2r3Ld 5d^ , 2nd^ 23L^ , 19, 

^ n 2 3L 0 ^ 2 24^ 

24n 



For the left-hand segment. 



,n,2,2Ld d^ 2nd^, 
^L'' ^ n 2 ■ 3L ^ ' 



for 0 < d < 

— — 2n 



^DtS(l) 



19 , ,n, 2 , 3 Ld - 5 d^ , 2 nd^ 23 L^, ^ L ^ ^ ^ l 

24 2 -^^L-'— 2 ^' IK ^ 

24 n “ 



1 ; 



for d > — 
— n 



(5) 



2 . Middle Line Segments 

We next develop the cumulative distribution functions 
for the minimum of the distances between a random trouble spot 
and n-2 carriers assigned to patrol stations i = 2,3, ..., n-1; 
stations which have at least one other patrol unit on either 
side of their patrol segment. 
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The middle line segment conditional cdf is identical 
for each of the "middle" patrol stations and is distinguished 
from that associated with either of the end patrol stations 
as follows: 

(1) end segments 

D = min(D, ,D») or D = min(D , ,D ) . 

12 n-1' n 

(2) middle line segments 

D = min for i = 2,3, ..., n-1 . 

L 2L 

A typical middle line segment is segment 2, ( ^ ' n^ * 

The assumptions for developing this portion of the model 
are : 

(1) T ~ Uniform C0,L), so T ^ with probability 

1 

n' 

(2) ~ Uniform (0/~) r 

(3) X- ~ Uniform (^/-%)/ 

^ n n 

(4) X, ~ Uniform (— ,— ) 

o n n 

(5) X^,X 2 ,X 2 / and T are independent random variables. 

a. Conditional Probability Density Function 
Consider a given t e 



Di = |t-X^| , 

D2 = |t-X2l, 

D = min{Dj^,D2) , 
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2L 

and the maximum value of D is 1. For 0 < d < t 

n — — 



n 



^D| T,S (2) ' P [D e (d,d+Ad) ] 



= P[X 2 £ (t-d,t-d-Ad) or X 2 e (t+d,t+d+Ad) ] = 2Adn/L. 



For t - — < d < — - - t , 
n — - n 



P [D e (d, d+Ad) ] = P[{X^^e (t-d,t-d-Ad) and ^2 > (t+d+Ad) } 



or {X 2 e (t+d, t+d+Ad) and X^ < (t-d-Ad)}] 



- (t+d+Ad)] • ^ ^ [t-d-Ad]-^ 

Li n Li Li Li 



2Adn rL 



- d - Ad] 



2Adn 



[1 - 



(d+Ad)^] 



3L 2L, 



Consider next a given te ^2 ~ I ” ^2 ' ' ^3 ~ 

D = min(D 2 ^D 2 ) and the maximuin value of D is t 



t-X3|, 
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For 



P [D £ (d,d+Ad) ] 



For 



P [D e (d,d+Ad) ] 



Thus it is seen 



^d|t,S(2) 



0 < d £ ^ - t , 



= P[X^ £ (t+d,t+d+Ad) or £ ( t-d , t-d-Ad) ] 



2Adn 

L • 



^ -t < d < t - 



P[{X 2 £ (t-d,t-d-Ad) and X^ > (t+d Ad)} or 



{X^ £ (t+d,t+d+Ad) and X 2 < (t-d -Ad) }] 



^ . [It- (t+d+Ad) ]^+ ^(t-d-Ad-t)t 

L n L L n L 



2Adn rL , . 

— - ■J - ''-il 



2Adn 



[1 - (d+Ad)n/L]. 



that : 



2n 
L ' 



2nn dn. 



0 < d < t- 5 



n — — n 



^ ,L 3L> 

'h '2ir> 



2n 

L 



2 n r T 

T 



0 <_ t 



- t<d<t-i 

n — n 



. ,3L 2L. 

t £ 

2n n 



( 6 ) 
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L 3L 

b. Unconditioning Over t e (— /-^) And 
. 2 t n 2n . 

t £ Given S(2) 

The conditional probability density function is 
graphed below, where values of fj^|,p(d|t) are given for 
various regions of the d-t plane. 




Figure 4 . Unconditional probability density function 



The 

is now found by 



conditional distribution of D, given S(2), 
taking the expectation of s(2) with 
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respect to the conditional distribution of T given, S( 2). 
Thus, unconditioning over T (given S(2)) has the following 
mathematical form: 



^D|S (2) 



2L/n 

^ ^D[T,S (2) ^T|S(2) 

L 



n 



For 0 _< d ^ , 



d+(L/n) 



^D| S ( 2 ) 



/ 



2 n I n X f n \ j f ■ ( ^ \ j 

^(1 -;f) (T-)dt + / 



2L 

n 



d' 'L 



L/n 



d+ii 

n 



2n ,n. 



L 'L' 



2L/n 



(2L/n) -d 



,n> 2 , , nd 



L L 



n> 2 , 2L 



= 2(^)"(l-i^) (d + ^-^) + 2(^)^(d^- d - d - ^ ) 



n n 



n 



n 



+ 2 ^ (1 - - — + d) 

L L n n 



,, ,n^ 2 ,L 2nd 
^ '■n “ L * 



for i < d < §, 



^D|S(2) 



(f-d) 



L 

n 



2L 

^(1 - (g)dt + / ^(1-nd) (g)dt 

dt(|) 
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^D| S (2) 



2 ( 2 .) 2 ( 1 - 1 ^) 

L L n n 



+ 2 ( 2.)^(1 -^)(^^- d -^) 
L L n n 



2(^)^(— - 4d + ^^^) . 

L n L 



Therefore in any middle line segment, 



T ,n^ 2 ,L 2nd , 



0 < d < 



L 

2n 



'D I S (2) 



(d) 



2(g)^(^ - 4d - 
ij n 



2nd^ . 





(7) 



c. The Conditional Cumulative Distribution Function 
The cdf of D conditioned by S(2) is calculated 
by integrating the pdf over the range (0,d) of D; 



^d| S (2) ^ ^ 

0 



Di S (2) 



(h) dh 



For 0 _< d , 



^DlS(2) 






2nh' 



-) dh 



^ ,n. 2 Ld 2nd . 
^ n " 3L ^ 
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and 



F ( ^ ^ 

'^D|S( 2 ) 6 






^D| S ( 2 ) 



2nh 



|. / 2(a,2,2i-.4h. , 

Li 

2ri 



•) dh 



5 ^ /Hv 2 r 2Ld. j 2 I 2nd. ^ ^ ^ i ^ ^ ^ ^ 

?*2(j;) [--2a + TIT ■ <— > ‘2H> ■" 2 <2H-> 



- (—) (i) ^] 

^ 3 L^ ^ 2 n^ ^ 



5 , .,,n,2r2Ld ^^2 , 2nd^ 7L^ , 

6 ^ 2(^) [--2d + _-^]. 

12n 



For the middle line segment, 



^D| S (2) 



,n^ 2 ,Ld 2nd ^ 

2‘l> <— ‘Tir’’ 



for 0 £d 



5 , ,n, 2 ,2Ld ^^2 , 2nd^ 7L^ , ^ L ^ ^ ^ 

6-'2(l) (--2d +^L TT)' 

l2n 



1; 



for d > — . 
n 



( 8 ) 



3 . Right-hand Segment 

Development of the conditional cumulative distribution 
function, of the minimum distance for the right-hand segment is 
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P |t-< 



identical in methodology and form to the results obtained 
in analyzing the left-hand segment in paragraph II.B.l above. 
We therefore omit the details. 

C. UNCONDITIONAL CDF FOR n SHIPS ON A LINE OF LENGTH L 

The cumulative distribution functions derived thus far 
are conditioned on S(i), i = 1,2, ..., n. Since T - Uniform 
(0,L) and each patrol segment is identically ^ units long, 
the probability of T occurring in any specific segment is , 

P[S.] = p[T £ = i. 

1 n n n n n 

Therefore the unconditioned cdf is derived from (5) and (8), 
as follows: 



n 



F^(d) 



= - I P 

n ^ 



D| S (i) 



i=l 



For 0 £ d ^ , 



Fo(d) 



2 3 3 

/ 1 \ r r / \ ^ ^ Ij^i 2 n » -i , '^ \ r / \ ^ / o \ / 2 n \ i \ 

(-)(2[(j;) + (n-2)[(j;) (2)(— 



2 2 3 3 3 

2nd nd 4n^d - 4n^d 



nL^ 



3nL' 



6L^nd - 3Lnd^ + 4n^d^ - 4n^d^ 
3L^ 
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and 



F (— ) 
D 



6 



12 



n = 



lOn - 1 
12n ■ 



For ^ < d < 



Fj^(d) 



(lOn-1) . ,l^ r-, r 2 ,3Ld 5d^ 2nd^ 23L^^ 

T2H 2 3ir~T~I^ 

2 4 n 

, , ,n^ 2 r2Ld -,j2 , 2nd^ 7L^ i 

+ (n- 2 ) ( 2 ) (:jr) [— 2 d +—^57 7 ]) 

L n 3L 3_2n^ 



(lOn-1) 4nd-2d 3nd^-4n^d^ 4n^d^-4n^d^ 5 7 

12n L ^2 ~ ^ 



22 2 22 32 43 333 3 

12n L d-6L nd+9n^d L-12n'^d^L+4n^d-^-4n''d'^+L''-nL'" 

or 3 

3L n 



Therefore the ciimulative distribution function of the 
minimum of the distances from T to n patrol units is: 



2 2 2 3 3 3 

6L nd-3Lnd +4n d -4n d 



3L' 



for 



0 < d 



< 



2n 



Fo(d) 



22 2 22 32 43 333 3 

12n L d-6L nd+9n d L-12n d L+4n d -4n d +L -nL 



3L^n 



for 





( 12 ) 
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III. APPROXIMATIONS TO ANALYTICAL RESULTS 



A. RATIONALE FOR APPROXIMATIONS 

The existence of a simple, readily applicable approxima- 
tion, especially one with well understood and widely known 
characteristics of behavior, might be useful to the military 
commander and planner. This thesis investigates two such 
approximations and evaluates the range of values for which 
each may be an appropriate approximation of (12). 



B. POISSON PROCESS 

The Poisson process as a mechanism for positioning n 
units on a line of length L is appealing because of the well 
developed and tractable body of theory regarding the exponen 
tial distribution of the distances between units. The limi- 
tations of the approximation will be discussed after the 
model has been completely developed. 

1. General Properties 

For the model let n be the average number of ships 
located on a line of length L, so A = ^ is the average den- 
sity of ships per unit length of L. Let be the distance 
between unit i and i+1, for all i; i= 1,2, ..., n-1. Then 

for a general Poisson process one has; 



P [N = n] 



/ . • \ ri All 
( AL) e 

n! 
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The interarrival "distances" (W^^) between adjacent ships are 
independent exponentially distributed random variables such 
that 



- Exp (1/A). 



In what follows only those cases where T occurs between 
two units will be considered (the interval (0,X^) will be 
discussed later) . Let the trouble spot T occur at random 
on the line of length L as illustrated below. 

K K H 

J / 

K— A12. H 




Let Ml represent the distance, from to t, M2 the distance 
from t to ^^en it is easily seen that Ml + M2 = 

the distance from the i^^ unit to the i+1^^ unit. The 
memoryless property of the exponential distribution implies 
that both Ml and M2 are distributed as exponential (1/A) 
random variables. Now with D = min (Ml, M2) , we have. 



Fj^(d) = PlMin(Ml,M2) £dj = P(.D^d) 



= 1 - P (D >d) = 1 - P CMl > d,M2 ^d) 



= 1 - P (Ml > d) P (M2 > d) 
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p 





and since Ml and M2 are identically distributed, 

F^{d) = 1 - [1 - F^(d)]^. 

Thus one has, 

Fj^(d) = 1 - [1 - (1 - 

= 1 - e"^^^ , X > 0 . (13) 

The model developed above for approximating the 
distribution of the minimum distance from a random trouble 
spot to the adjacent patrol units implies the cdf. 

2nd 

Fjj(d) = 1-e 0<d (14) 

which is an exponential 2X cdf . 

2 . Limitations To The Poisson Approximation 

The Poisson approximation cannot guarantee that n 
units will be positioned within the (0,L) limit of the line 
to be patrolled. For example there exists a positive proba- 
bility, however small,- that no units will occur within 
(0,L) . 

The model as developed does not consider the specific 
distributions involved when T occurs to the left of the 
location of the first unit. The effect of this simplification 
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is to increase the size of d for an arbitrary probability 
p which results in a weakening of the boundary condition 
developed in III.D. 

The Poisson approximation does not preclude three 
or more units being very close to each other. 

C. FIXED STATION MODEL 

A second approximation, the fixed station model, is nov; 
considered. This model assumes each patrol unit maintains 
a fixed station located at the midpoint of the assigned 
patrol line segment = (2i-l)L/2n. 



o 



X, 




X4 

-1 ^ 

3L 

H- 



L 



1 . Model Development 

Since the distance between adjacent units is constant 

(^) and the distance from any unit to its patrol line segment 

boundaries is L/2n, the model developed for an arbitrary 

patrol segment will be representative of all segments. 

Consider segment 1. Given S(l), T is distributed 

uniformly on (0^^) ^rid the model being developed will specify 

the conditional cdf of D = min{-:^^ — T,T--^}. That is, 

Zn Zn 
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= P (D < d| S (1) ) = 



P(D < d|T <^)P(T;^^) 



+ P (D < d I T > 4 ^) P (T > 

— ' 2n — 2n 






P(T> A-d|T<^)i + p(T<d + ^|T>^)i 



. Ji L. 

In /rr, L j|_ L,, , Ir 2X1 2X1 , 

2 [1 - P (T £ 2n “ ^ 2 X 1 ^ ^ 2^ L/2xi ^ 



^ !<¥• 



In _ /L-2xid^ ,2n^ 1 ^ nd 
2 ^ 2n ^ ^ L ^ ^ L 



1 L - 2xid xid 

2 2L L 



2xid 



for 0 < d < 

— — 2x1 



(15) 



Expressioxi (15) is coxiditioxied by the occurrexice of 
S(l). However the coxiditional cdf givexi S(i) for each i; 
i = 1,2, . . . , XI ; of D = min (D^,D^_i_^) is identical to (15) 

and the summation of n conditional cdf's multiplied by the 
conditioning probability 1/n results in: 



P(D < d) 



2nd 
~ ' 



for 



0 < d 




(16) 



which has a median value 




(17) 
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2 . 



Limitations Of The Fixed Station Model 



The obvious limitation of the model is the loss of 
realism in constraining a mobile unit to a small geographically 
fixed station. As a consequence of this constraint the 
maximum value of D = min (D . ,D . ) is for the fixed station 

model . 

D. COMPARISON OF ANALYTIC AND APPROXIMATE MODELS 

The easiest method of contrasting the two approximations 
with the original analytic model is through a tabulation of 
their individual cumulative distribution functions for 
specific values of n, L, and D. In Table 1 d is tabulated 
as a percentage of the ratio This procedure allows one 

to readily compare the various models at differing L and n 
values . 

For the L = 30000, n = 6 case the distribution function 
for the analytic model is bounded below by that of the Poisson 
process and above by that of the fixed station model. 

Figure (5) is a graphic presentation of the "fit" of 
the 2 bounding models when compared to the analytic model 
with parameters L = 30000 and n = 6. The graph depicts the 
residuals obtained from differencing the analytic and 
approximation models at various values of d. Data points 
are from Table 1. 

For the case under consideration the Poisson model median 
value of D (0.3465) does not approximate the analytic model 
median (0.266) as closely as does that of the fixed station 
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EVALUATION OF ONE DIMENSIONAL ANALYTIC AND APPROXIMATING MODELS 
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Table 



FIGURE 5. Plot of Residuals of Approximating Models 
Versus Analytic (n = 6, L = 30000) model 
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model (0.25). Except for values in the neighborhood of 
d = .5 the Poisson process model residuals are of greater 
magnitude than are the fixed station residuals. 

The Poisson model gives reasonable (if somewhat conserva 
tive) estimates of the analytic model results within the 
range 0.0 _< d j< .3. 

For this case, the fixed station approximation model 
not only bounds the analytic model, from above, it also 
appears to be a reasonable approximation of the analytic 
model over most of the range of de (0.0, 1.0). The fixed 
station median value of 0.25 is very close to the analytic 
model median value of 0.266. Using a criterion that a 
residual value greater than 0.1 indicates that the approxima 
tion model is inappropriate for that value of d, it is seen 
that the fixed station model is appropriate for all d except 

0.4 < d < 0.6. This result, while developed for the case 
n = 6 and L = 30,000 will hold in general except in those 
unusual cases where n > L. 

The results tabulated in Table 1 suggest that there 
exists a natural hierarchy or ordering of the three models. 
If one fixes the P(D^d) = .5 and then ranks the models 
according to the relative magnitudes of d required by each 
model to obtain this probability the result is as follows: 



Model 



Median 



1. Poisson 



0.3465 



2. Analytic 



0.266 



3. Fixed station 



0.25 
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Arranged in this manner the models reflect the following 
"costs" : 

1. Increasing level of control by higher operational 
authority 

2. Stochastically decreasing distances to random 
targets 

3. Decreasing randomness of own units motion 

4. In a hostile environment an increasing risk of an 
adversary accurately determining the positions of 
the patrol units 

5. Decreasing probability of adjacent units being in 
close proximity (overlapping of coverage) which is 
paralleled by decreasing problems with mutual inter- 
ference and contact identification. 

The geometry of the residuals plotted in Figure 5 indi- 
cates that the average of the two boundary cdf ' s might form 
a more reasonable estimate of the analytic model than either 
of the boundary models does when considered individually. 

The boundary model data in Table 1 was averaged, sub- 
tracted from the analytic (L = 30000, n = 6) model and the 
residuals plotted on Figure 5 as the average model residuals. 
The resulting residuals are obviously a very reasonable 
conservative approximation and lower bound for the analytic 
model, particularly for those de (0,.65). 
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IV. APPROXIMATIONS FOR TWO-DIMENSIONAL PATROL AREAS 



The previous chapters have developed a one -dimensional 
patrol model and two approximations which form upper and 
lower bounds for the analytic model. It seems a reasonable 
supposition that when the problem scenario is expanded to 
include units patrolling a two-dimensional area that a two- 
dimensional Poisson field would form a lower bound of the 
two-dimensional analytic model and that a two-dimensional 
lattice of fixed stations could reasonably be postulated as 
forming an upper bound on the two-dimensional analytic model. 

The algebra involved in developing the two-dimensional 
analytic model for n randomly patrolling units is quite 
involved and time constraints precluded its development in 
this thesis. However the two-dimensional approximations are 
tractable and are presented here. 

A. TWO-DIMENSIONAL POISSON FIELD 

1 . General Theory For A Poisson Field 

Consider an array of points distributed over a 
plane P. For each subset B of P let N(B) be the number of 
points within B. The array is said to be distributed in 
accordance with a Poisson process with density ? if the 
following assumptions are met: 

(a) the number of points in non-overlapping areas 
are independent random variables 

(b) for any subset B of finite area, N(B) is Poisson 
distributed with mean ^b, where b is the area 
of B. [Ref. 1] 
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The dimensional analysis of ;; is — 



length, squared' 

Thus if N(b) represents the number of points in a region 

with area b; 



P(N(b) = n) = 



(gb)^ e 
n! 



( 20 ) 



The probability that there are no points within a radius r 
of a given location is determined by 

2 

P (R > r) = P(no points within a region of area irr ) 



= e 



-TTr C . 



( 21 ) 



2 . Poisson Field Model 

a. Assumptions and Symbology 

The following assumptions are made: 

(1) The n units are distributed stochastically 
by a Poisson process with rate 

(2) The trouble spot T is equally likely to 
occur anywhere within the patrolled area. 

(3) The area to be patrolled is L units in 
length and W units in width (L > W) . 

b. Development 

In accordance with the theory for a Poisson 
field, C = ^ and following expression (21) , the Poisson 
field model is. 
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P (D < d) 



= 1 



-Trd 



2 



e 






= 1 - e 



-TTd n/LW 



( 22 ) 



c. Limitations to the Poisson Field Model 

The Poisson field model cannot guarantee that 
n units will be positioned within the patrol area. Nor can 
the model preclude five or more units from being assigned 
positions which are in very close proximity to each other. 

B. FIXED LATTICE OF POINTS 

The cdf of the model which formed the upper bound of the 
one-dimensional analytic model was a series of patrol 
stations fixed in the middle of each patrol line segment. 

A reasonable two-dimensional approximation would seem to be 
an array of patrol stations fixed in the centers of their 
respective patrol areas. It is hypothesized that this two- 
dimensional array is an upper bound for the two-dimensional 
analytic model. 

1 . Assiimptions 

Assume n points are located within an L by W (L > W) 
area as shown in Figure 6 below. 

Assume that the location of T is uniformly distri- 
buted over the L by W rectangle. The n ships are to maintain 
station at the center of their respective L//n by W//n 
areas . 
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Figure 6. Fixed lattice of points 
2 . Development 

Since each of the n ships is centered in an iden' 
tically shaped area the model developed for any one unit 
will hold true for all units . The maximum value of D is 
the distance from the center to the corner of a unit's 
patrol area as shown in Figure 7 . 

a. Conditional Cumulative Distribution Function 
Define event S(i) as the event T e (region i) 
with P[S(i)] = 1/n, and 

^DjS (i) ^ (D < d|S (i) ) . 
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Patrol Area of Fixed Lattice Model 



- 4 . 




length d d, the radius of the expanding circle, 

9 = 2 cos"^ ^ = 2 cos"^ 

arc length s = d0 , arc length t = dcj) 

A3 = arsa of the segment under s but outside the rectangle 

= i(ds-ga) = d^ cos'^(-iL_, - (^2 ,W^)l/2, W 

A - 2d/n 2/n 

area of segment under t but outside the rectangle 

= (dt - he) = d^ cos — ) - — ^ ^1/2 

2d/E 2/n 

Figure 7 
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2 /n 
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1 ^^ 




Refer to Figure (7) and consider the case of 
the inscribed circle of radius d such that, 0 < d < W/2/n, 



P(D<d|S(i)) = 



n 



nird 

LW 



Consider the case of W/2/n £ d £ L/2/n where 
the expanding circle has intersected the longer sides of 
the patrol rectangle. 



I S (i) 



- (2A ) ( ^ ) 

LW ^ s^ ^LVl/n’ 



As the circle radius continues to increase, consider 



L/2/K < d < « 



for which 



^D| S (i) 



M|-- (2A^ + 2At»<Lvk' 



Thus the conditional cdf is , 



^D| S (i) 



mrd 

LW 

nird'^ 

LW 

mrd^ 

LW 



2nA 

C 

“lw" 

2nA 

C 

~W 



for 0 < d < 



W 



2/n 



2nA^ 

Tw" 



for 



for 



w 

2/n 
L 



< d < 



2/n 

,2 



2 /n 



< a < ^ ^)i/2 

— — Un 4n 
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b. Cumulative Distribution Function 



To uncondition over the event S(i) one needs to 
take the following summation, 



n 



FoCd) 



n ^ ^D[S(i) 
i=l 



^dI S ( i) ' 



since the n conditional cdf's are identical. Therefore it 
is seen that: 



Fo<d) = 



nud 
LW ' 

nird^ - 2nA 



LW 

nird^ - 2nA 

s 



for 0 < d < 



W 

2/n 



2nA. 



LW 



for 



for 



w 



2/n 

L 

2/n 



< d < 



2/n 

,2 



^ /W ^ L^, 1/2 

< d < (t— + -j—) 

— — '4n 4n' 



(23) 



C. EVALUATION OF THE APPROPRIATENESS OF TWO-DIMENSIONAL 

APPROXIMATING MODELS 

Expressions (22) and (23) form hypothetical bounds to 
an undeveloped area patrol model. If the hypothesis is 
correct it would seem reasonable that the two-dimensional 
bounding functions should behave in a manner similar to that 
of the one-dimensional bounding functions. That is one 
would not only expect the fixed array model to yield sto- 
chastically larger results than the Poisson field model but 
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one would also expect to see very close agreement between 
the two models for small values of d. Table 2 which tabu- 
lates the evaluation of the two models for parameters n = 25, 

L = 30000 and W = 20000, indicates that for this case the 
two-dimensional models are indeed well behaved. 

In the one-dimensional case it was found that the average 
of the boundary models was itself a very reasonably estimator 
of the analytic model. If the hypothesis regarding the two- 
dimensional boundary models is correct then a similar averaging 
process could be considered a better approximation of the 
analytic area patrol model than either boundary model is 
when considered separately. The average of the two-dimensional 
boundary models is also tabulated in Table 2. 



Model d values* (in units of length) 





0 


500 


1000 


1500 


2000 


2500 


3000 


3500 


3605 


Poisson field 


1 


.032 


.123 


.255 


.408 


.559 


.692 


.798 


.818 


fixed point 
array 


1 


.033 


.131 


.295 


.524 


.733 


.920 


.998 


1.0 


average 


1 


.032 


.127 


.275 


.466 


.646 


.806 


.898 


.909 



* for parameters n = 25, L = 30,000, W = 20,000. 



Table 2. Evaluation of Two-dimensional 
approximating models 
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V. ANALYTICAL MODEL WITH RADIUS OF INFLUENCE 



The initial analytic model, equation (12) , was developed 
under the assumption that an aircraft carrier be considered 
a point on its patrol line segment. This restriction effec- 
tively eliminated from the previous analysis the known ability 
of a carrier's airwing to exercise considerable influence 
over the environment adjacent to the carrier. We now modify 
the model to allow an examination of the contribution to 
P(D_<d) made by an airwing controlling a radius r on either 
side of the carrier's position. 

An additional benefit of relaxing the zero radius of 
influence assumption may be a quantitative model of how a 
units' probability of being within ,a specified distance of 
a random target will vary as the radius of influence is 
changed . 

A. ASSUMPTIONS^ 

All the assximptions made for the initial analytic model 
remain in effect with the exception of assumption II. A. (5) 
which constrained a carrier to being considered a point on the 
patrol line segment. In this modification each carrier will 
be considered to have a radius of influence, r, on either side 
of the carrier's position. If the trouble spot T occurs within 
± r of the carrier's position the distance from the carrier 
to the trouble spot is defined as zero. Distance to the 
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trouble spot, T, will be measured to the edge of a carrier's 
radius of influence instead of to the carrier's actual 
position as was done in Chapter II. 

There are no restrictions regarding the movement of a 
carrier along its patrol line segment. This means that the 
second carrier would not be restricted from periodically 
exerting influence over a portion of the first r miles of the 
third carrier's line segment (see in Figure (8)). 

B . MODEL DEVELOPMENT 

1. Cumulative Distribution Function 
In the initial analytical model, 

D = min(Dj^,D 2 , . . . ,D^) where De * 

In the current modification 

H = min(Hj^,H 2 / . . . ,H^) where H z (0, "r) 

and since D^ = + r for i = 1,2, ..., n then 

D-r for r < D < — 

H = ^ 

0 for 0 < D < r. 



Now for h > 0 , 
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t h 

the position of the i carrier 

the radius of influence on either side 
of a carrier's position 

length of the line to be patrolled 



n 



the number of carriers patrolling 
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FIGURE 8. Schematic and Symbology for the Radius of 
Influence Model 
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Fy(h) = P[H£h] = P [H £h|D ^r]P (D ^r] 

+ P[H£h|D£r]P(D£r] 

= P [D-r £ h|D ^r]P [D ^r] + P [0 _< h | D £ r] P [D £ r] 

_ P [D-r < h,D £ r] P [D £ r] P [ 0 £ h , D £ r ] P [D £ r] 
P[D£r] P[D£rj 

and since [0 £ h] is certain, the above reduces to 
P[r£D£h+r] + P[D£r] = Fp(h+r) - 

= Fj^(h+r) ; 0 £ h £ | - r. 



(24) 



The radius of influence model is 

Fjj(h) = Fjj(h+r) for 0 £ h £ ^ - r (25) 

which is the initial analytic model [expression (12) ] 
evaluated at D = h + r. 

2 . Marginal Rate of Return With Respect to Radius of 
Influence Model 

When the radius of influence model is differentiated 
with respect to r the result is the equivalent to evaluating 
the pdf of equation (12) at h+r: 

f^(h) = f^(h + r) (26) 
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At r = 0 this equation is the true pdf associated 
with expression (12) ; at values of r > 0 this function gives 
the expected marginal retvurn in terms of P(H£h) for a unit 
change in the radius of influence. 

C. MODEL EVALUATION' 

The following plots were developed to assist in visualizing 
two specific realizations of the cumulative distribution 
function and marginal rate function for the radius of 
influence model. 



n 



Type of plot 



Figure 



10 


10 


cdf 
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10 


10 


marginal rate 


10 


30000 


6 


cdf 


11 


30000 


6 


marginal rate 


12 



The parameters L = 30,000 miles and n = 6 ships were 

chosen as being representative of the imaginary but roughly 

realistic situation described in the Introduction. 

The plot in Figure 11 appears to be quite linear for 

all values of r up to and in some cases beyond the r = 0.0 

median value of 1330 miles. It is also of interest to note 

that for large ratios of il , that for a given h a constant 

n 

increase in r results in a near constant increase in P(H_<h). 
This characteristic is not observed in the L/n = 1 case. 

By the time the distance from the target has reached 
3500 miles the n = 6, L = 30000 cdf has exceeded 0.95. In 
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other words, for all values of r ^ 0 P(H^3500) > .95. 

While this may be a counter-intuitive result it is not 
unreasonable if one remembers that a value of h = 3500 miles 
encompasses a 7000 + 2r mile portion of the total line. 

D. EMPLOYING THE RADIUS OF INFLUENCE MODEL 

We consider some examples which illustrate use of the 
model. The operational commander may desire to know the 
minimum number of units necessary to patrol a 30000 nautical 
mile coast such that the nearest unit has at least a 60% 
chance of being within a 2-j day transit of a random spot 
along the coast. Further he may desire to know if there 
is any substantial advantage to loading the carriers with 
Type A aircraft (180 nautical mile radius) as opposed to 
loading with Type B (300 nautical mile radius). The units 
steam at 25 kts . 

Solve equation (25) for n at h = 1500, r = 180, L = 30000 
and P(H_<h) ^ .6. The result is n ~ 6. Equation (25) , 
evaluated at n = 6 , L = 30000, h = 1500, r = 180 results 
in 

F„(h) = 0.611 > 0.6, 

rl — 



as required. 

The marginal rate of return function (26) , with the 
above conditions gives 
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f„(h) = 0.0003, 

n 

which indicates one could expect a 0.0003 increase in 
F„(h) for each unit increase in r. The Type B aircraft's 

n 

300 nautical mile radius will increase F„(h) by approximately 

n 

(300 - 180) (0.0003) = 0.036 over the value of 0.611. (The 
actual increase by solving expression (25) for both R = 180 
and R = 300 is 0.03959.) 

As a result an operational planner would inform the 
commander that it would require a minimum of six carriers 
to patrol the line. Additionally, based on an MOE of 
P(D£d), there is no substantial advantage to be gained by 
loading a particular type aircraft on the six carriers. 
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CDF FOR FIXED VALUES OF' 
RADIUS OF INFLUENCE 
L= 10 N= 10 
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FIGURE 9 
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FIGURE 10 
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CDF FOR FIXED VfiLUES OF 
RfiDIUS OF INFLUENCE 
L= 30000 N= 6 
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FIGURE 11 
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FIGURE 12 
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VI. EXTENSIONS OF THE ANALYTIC MODELS 



During the course of the development of this thesis 
several potentially interesting extensions were discussed. 
Some of these ideas are presented below as possible areas 
of future investigation. 

A. MULT I -DIMENSIONAL PATROL AREAS 

The assumption limiting a carrier to patrolling along 
a line is unnecessarily restrictive. A natural extension 
would be to expand the patrol line to a patrol area of 
length L and width W such that each carrier would be assigned 
a L/n X W patrol area as illustrated below: 



T 

w 

1 



T, the trouble spot, would be restricted to either occurring 
along the (_0,L) line, or perhaps to occurring anywhere within 
the (L,W) area. 

An extension which moves even farther from the patrol 
line model, but which naturally follows the previous exten- 
sion is to consider a grid imposed on an open ocean area with 
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one carrier assigned each rectangular patrol area. T would 
occur randomly anywhere within the patrol area boundaries. 




W 



Figiire 13. n carriers randomly patrolling a 2 dimensional area. 

The radius of influence of each carrier could also be 
incorporated into the two-dimensional models . 

B . GAME THEORY 

The stipulation that the location of the trouble spot 
is a random variable with a particular probability distribu- 
tion may be realistic in most situations involving a naval 
presence mission. Where an adversary relationship can be 
postulated this stochastic assumption as well as the assximp- 
tion of independence between the location of T and the 
positions of the carrier units (X^'s) may be untenable. If 
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the adversary is capable of some degree of knowledge of the 
carrier's positions, through radio direction finding, 
shadowing or satellite reconnaissance, he is capable of 
choosing T so as to maximize his return (perhaps by maxi- 
mizing the expected distance to the patrolling force) . In 
this case the problem takes on a game theoretic aspect, and 
results from game theory might be applicable. 

If, through careful analysis of an adversary's intelli- 
gence gathering capabilities, force levels and historic 
political intentions, one is able to determine an a priori 
distribution of the location of T, then it is possible to 
use the probability integral transform to transform the a 
priori distribution into a Uniform (0,1) distribution. This 
might enable one to utilize the models developed herein 
and then transform the results back to the original a 
priori space. 

C. ECONOMIC COST OF VARYING THE RADIUS OF INFLUENCE 

There are economic costs involved in varying the radius 
of influence of a carrier. Through the use of the radius 
of influence model developed herein and basic cost considera- 
tions, an analysis of the cost tradeoffs involved in achieving 
a particular improvement in P(D£d) could be developed. This 
development could include economic tradeoffs in terms of 
various combinations of: 

1. Increasing a carrier's radius of influence 

2. Increasing the transit speed of the carrier 
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3. Increasing the number of carriers 

4. Increasing .intelligence concerning the location 
of T. 
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